On the nature of the magnetic ground-state wave function of V2O3. 
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After a brief historical introduction, we dwell on two recent experiments in the low-temperature, 
monoclinic phase of V2O3: K-edge resonant x-ray scattering and non-reciprocal linear dichroism, 
whose interpretations are in conflict, as they require incompatible magnetic space groups. Such 
a conflict is critically reviewed, in the light of the present literature, and new experimental tests 
are suggested, in order to determine unambiguously the magnetic group. We then focus on the 
correlated, non-local nature of the ground-state wave function, that is at the basis of some drawbacks 
of the LDA+U approach: we singled out the physical mechanism that makes LDA+U unreliable, and 
indicate the way out for a possible remedy. Finally we explain, by means of a symmetry argument 
related to the molecular wave function, why the magnetic moment lies in the glide plane, even in 
the absence of any local symmetry at vanadium sites. 
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I. INTRODUCTION 

In the last 35 years V2O3 has been the subject of 
wide experimental and theoretical investigations, af- 
ter its identification as the prototype of Mott-Hubbard 
systems*. At the beginning of the seventies an inten- 
sive set of measurements (resistivity 2 , specific heali*, 
magnetic susceptibility*, x-rays*, and polarized neutron 
diffraction 6 -) was performed in order to clarify its crystal 
and magnetic structure as well as its transport properties. 
This marked a break-point with all previous theories** 
that described the "metal-semiconductor" phase transi- 
tion of V2O3. Actually, two metal- insulator transitions 
were found: one from a paramagnetic metallic (PM) to an 
antiferromagnetic insulating (AFI) phase, below ~ 150 
K, associated with a corundum-to-monoclinic crystal dis- 
tortion; the other from the PM to a paramagnetic insu- 
lating (PI) phase at higher temperatures (above ~ 500 
K): among transition metal oxides it is the only known 
example of PM-PI transition*, with no structural defor- 
mations and no magnetic degrees of freedom involved. 
The PI phase can be also obtained by means of a small 
(1 -j- 2%) Cr-doping: in this case Cr ions act as a neg- 
ative pressure, thus enlarging the average cation-cation 
distance and inducing a Mott-insulator transition. As 
described ini*, all transitions are due to the interplay 
between band formation and electron Coulomb correla- 
tion. 

In 1978 Castellani, Natoli and Ranninger (CNR) at- 
tempted, in a series of papers^*, a realistic description 
of the complex magnetic properties and phase diagram 
of V2O3. CNR focused on the peculiar magnetic struc- 
ture observed in the AFI phase, that breaks the high- 
temperature trigonal symmetry and can not be explained 
in terms of a single-band Hubbard model. They real- 
ized that the introduction of the extra degree of free- 
dom represented by t2 S -vanadium orbitals was necessary 
to describe correctly the spin structure. The driving 
mechanism was the ordered pattern of these t2 9 -orbitals 
throughout the whole crystal, the so-called "orbital or- 
dering" . Twenty years later, the advent of modern third 



generation synchrotron radiation sources stimulated an 
interesting prediction*- to look for the expected orbital 
ordering by means of resonant x-ray scattering (RXS). 
Soon after, its experimental detection was claimed to be 
founcU*. The subsequent theoretical and experimental 
debate led to the abandonment of the CNR model, that 
was based on a spin Sy = 1/2 on each vanadium ion 
and thus in contradiction with the experimental result of 
Sy = l*i*i*, in favour of a picture where the two near- 
est V-ions are linked together in a stable molecule with 
spin Sm = a*** 6 *!*. This picture has been recently crit- 
icized by Elfimov et alM.. Yet, the LDA and LDA+U 
schemes at the basis of this work are unable to describe 
the physics of V2O3, due to a fundamental weakness of lo- 
cal density approximations to mimic the behaviour of the 
correlation potential, as we shall see in Section 4. Thus, I 
believe that, especially in the version of Tanakai* that in- 
cludes the spin-orbit interaction, the "molecular" model 
catches the main energy scales and correctly describes 
many relevant ground state properties: spin magnitude 
and direction, softening of the C44 phonon mode, spin 
to orbital moment ratio and relative orbital anisotropics 
(dig to e g ratios) in the three PM, PI, and AFI phases. 



The paper is divided as follows: in Section 2 the im- 
plications of the RXS experiments**!* about the nature 
of V2O3 magnetic space group are introduced. Section 
3 is devoted to a different interpretation of the "non- 
reciprocal" linear dichroism experiment in Ref. 20 , in or- 
der to reconcile it with the magnetic space group as de- 
termined by RXS. Finally, in section 4, the breakdown 
of both local density approximation and LDA+U is illus- 
trated, due to the non-local correlation potential gener- 
ated by the molecular ground-state wave function. The 
molecular ground state is also at the basis of a simple 
explanation, based on symmetry arguments only, of the 
reason why the magnetic moment should lie in the glide- 
plane, in spite of the absence of a local symmetry element 
at vanadium sites. 



II. INTERPRETATION OF RESONANT X-RAY 
SCATTERING EXPERIMENTS. 

The RXS experiments^ 3 -^ were predicted to occur on 
the basis of CNR model, and had of course to be rein- 
terpreted once the failure of such a model had been 
recognized. Such a task was undertaken by Tanakaii 
and Lovesey and collaborators^, who, independently, 
identified the RXS signal as a signature of the mag- 
netic distribution and not of the orbital ordering of 
vanadium electrons. Their interpretations are based on 
the commonly accepted magnetic space group (MSG), 
P2/a + T{E\to}P2/a, derived from x-ray^ and polar- 
ized neutron^ experiments, that contains explicitly the 
time-reversal symmetry T. Here to is the body-centered 
translation and P2/a the monoclinic group containing 
the identity E, the inversion /, the two-fold rotation 
about the monoclinic b m -axis C 2 b and the reflection rhb 
in the plane perpendicular to this axis. Of course, the 
appropriate translation is associated to each of these op- 
erators, as detailed below: 
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FIG. 1: Four vanadium ions lying in the a m -c m plane are 
shown. They are all ferromagnetically-coupled and related 
one another by the four symmetry operations of the P2/a- 
group. The dashed ellipse marks the vertical molecule. 



1) E, I 

2) C 2 b, rh-b 

3) f , ft 

4) fC 2b , fm b 



No translation 

\{o-m + b m + Cm) = to 



A = (e ia + Ie- ia )(l - Tm b )(l - T)/i 



(2) 



The origin of the axes has been chosen in the inversion 
center shown in Fig. ^ Referring to the figure labels, Vi 
and V2 are connected by the C2f,-rotation through their 
midpoint, V2 and V3 are related by inversion and Vi 
and V3 by the glide-plane operation rhb, within the plane 
of the figure. The unit cell contains eight V-ions: the 
other four are related to these by the time-reversal op- 
eration and the body-centered translation. Introducing 
the atomic scattering factor (ASF) as: 



the 

^8 



^Ku-(E n -E g )-iT n 
RXS amplitude can be written as: A 



(1) 



e l ®' Ri fj where the sum is extended over the 8 
equivalent vanadium ions of the magnetic unit cell 
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positions Rj, whose ASF is fj. Here ip g and tp n are the 
ground state and the photoexcited state, whose energies 
are E g and E n , respectively; the operator O represents 
the matter-radiation interaction, Huj is the photon en- 
ergy, r„ the inverse lifetime of the excited states and 
Q = (h,k,l) is the exchanged Bragg-vector. Reflections 
with h + k + l =odd and h =oddi2*iii are those of interest. 
Taking into account of the tensorial character of fji^LS^ 
and the properties of the MSG P2/a + f{E\t }P2/a, it 
is possible to express all the ASF in terms of just one of 
them, through the appropriate symmetry operations 23 . 
We get: 



where a = 2n(hu + kv + lw), and (u,v,w) are the 
fractional coordinates of vanadium ions in the unit of 
monoclinic axes. Equation Q has two important con- 
sequences: first, it shows that all reflections of kind 
h + k + l =odd and h =odd have a magnetic origin, 
in such a way that the factor (1 — T) be non-zero. In sec- 
ond place, it indicates that they are sensitive to magnetic 
multipoles of higher order than the magnetic moment, 
because of the factor (1 — Trhb). This latter observation 
stems from the fact that the magnetic moment lies on 
the glide plane rhb, thus being an eigenvector for Trhb, 
with eigenvalue +1. Such a picture clearly shows that 
the presence of the time-reversal symmetry plays a cru- 
cial role in the determination of the RXS origin. 

Yet, little before the publication o£i£ andSi, Goulon 
et al~ had interpreted the results of their linear dichro- 
ism experiment in the monoclinic phase of V2O3 as an 
evidence of the non-reciprocal gyrotropic effect, whose 
implications are in striking contradictions with the previ- 
ous conclusions. In fact, a necessary condition to detect 
the non-reciprocal gyrotropy tensor is that the system 
be magnetoelectric, implying, in turn, that neither time- 
reversal nor inversion symmetry can separately belong 
to the MSG. If the linear dichroism experiment had cor- 
rectly been interpreted, this would have implied a reduc- 
tion of P2/a+T{E\to}P2/a to a subgroup not containing 
f and I, like (E,fi,fC 2b ,fh b ) or {E,fi,C 2 b,Trh b ). 

Thus, one might wonder whether the RXS experiment 
could be explained with similar arguments as above even 
in the framework of the above magnetoelectric space 
groups, in order to be consistent with non-reciprocal ef- 
fects. My answer is negative, and is motivated by the 
results of Keism^^^^, that I summarize here. 
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For practice, we shall distinguish two physical mecha- 
nisms that can lead to a MSG reduction. The first is of 
"charge" origin, i.e., due to any kind of orbital ordering 
(not necessarily that of CNR); the second concerns the 
presence of an out-of-plane component of the magnetic 
moment directed along the b m axis. The first mechanism 
has been definitively ruled out by a numerical simulation 
of the RXS experiment, based on the finite difference 
method^ 7 -. In fact, in this case, due to the breakdown of 
the time-reversal symmetry in Eq. (J2J, for h + k+l =odd, 
h =odd reflections a signal would appear at 4p-energies 
much bigger than the one at the pre K edge, whereas 
no such a signal is detected experimentally. The second 
possibility has been long debated in the literatur e^ 2 ?' 2 ? , 
and yet is still very controversial. The importance of an 
out-of-plane component of the magnetic moment is re- 
lated to the fact that it could break the glide-plane sym- 
metry. Three independent polarised neutron scattering 
experimentsS* 2 ^ 2 ^ were not able either to discard or ac- 
cept this hypotesis. It is interesting to note here that in a 
crystal field approach, due to the absence of symmetry at 
vanadium sites, the magnetic moment would be allowed 
to have any direction. Nonetheless, this picture turns out 
to be too simple: the symmetry of the molecular wave- 
function is such as to force the moment in the glide plane 
(see Section 4). In any case, an out-of-plane component 
would again imply the presence of a signal at the (lll) m 
reflection in the energy region of the conduction p-band, 
contrary to what experimentally detected. In fact, as 
described in Refsii§iiL2i, the extinction rule that forbids 
this signal comes from the time-reversed glide-plane sym- 
metry, due to the factor (1 — Trhj,). This latter is active 
only when the magnetic moment lies within the mirror 
plane. For this reason, we can conclude that the absence 
of any features at 4p-energies oih + k + l =odd, h =odd 
reflections is an indication that the correct MSG of V2O3 
is P2/a + f{E\t }P2/a. 

Finally, I would like to underline a relatively minor, but 
still important difference betweenAI and2i. In Tanaka's 
paper the RXS signal is a consequence of the dipole- 
quadrupole (E1-E2) interference, thus sensitive to mag- 
netic, parity-odd physical quantities (like a toroidal mo- 
ment or a magnetic quadrupole - which of them is left 
unspecified). On the contrary Lovesey and collaborators 
were able to fit the experimental results of Refsii&iS in 
the different case of purely quadrupolar reflections (E2- 
E2), where the extracted information concerns the mag- 
netic octupole of the system. These two views seem con- 
flicting each other, but probably "in medio stat Veritas". 
In fact, a third more recent result on this subject^, based 
on ab-initio calculations, shows that both channels can 
contribute to the global signal: by means of a relativistic 
extension of the multiple scattering theory, implemented 
in the FDMNES package 27 , the authors of Rcf. 26 have 
been able to find a fair agreement in energy and az- 
imuthal scans with the experimental resultsi&i^. Both 
E1-E2 and E2-E2 channels are present in the global in- 
tensity and their relative contribution depends basically 



on the energy and on the azimuthal angle at which the 
reflection is measured. Yet, the single-particle character 
of these calculations still leaves open the door to future 
improvements in the description of the pre K edge RXS. 

III. INTERPRETATION OF THE LINEAR 
DICHROISM EXPERIMENT. 

As sketched above and already discussed in the 
literature 2 ^, the interpretations of RXS 17-21 and non- 
reciprocal dichroismSSi in their original forms require in- 
compatible MSG. Nonetheless, I believe that there is a 
way to explain both experiments within a unified theo- 
retical framework, which consists in a different interpre- 
tation of the linear dichroism experiment, as suggested 
in^. 

This alternative view is based on the possibility that 
the measured linear dichroism 20 . is not a gyrotropic sig- 
nal, but is related to the electric quadrupole, and thus 
of pure dipole-dipole (El-El) origin. The difference be- 
tween the physical operators associated with these two 
kinds of dichroism was first derived in Ref , 30 , where their 
formal expression was given. In a reference frame where z 
is the direction of propagation of the x-ray beams, and x 
and y, respectively, those of the associated polarizations, 
the pure El-El dichroism is described by the electric 
quadrupole operator L 2 — L 2 , while the non-reciprocal 
(E1-E2) signal is determined by the expectation value 
of (i 2 — Ly)n z — c.c Here L is the angular momen- 
tum operator and f2 is the toroidal moment2LS. The 
ab-initio simulations^, based on either the finite differ- 
ence method or the multiple scattering theory, clearly 
demonstrated the presence of a dichroic signal of El-El 
origin, the size of which is some percent of the absorp- 
tion signal. This calculated intensity has the same shape 
and order of magnitude as the one measured by Goulon 
et alm^. Indeed, apart from the ab-initio calculation re- 
sults, a dichroic signal of El-El origin and of this order of 
magnitude had to be expected, as the low-temperature 
distortion of the corundum hexagonal planes along or- 
thogonal directions, where the dichroism is measured, 
can be up to 4%. These simple considerations explain 
the findings of Refi^ also at a qualitative level. 

The reason why Goulon et al*& have interpreted the 
dichroic signal as "non-reciprocal" is mainly due to the 
fact that they found an inverted signal when their exper- 
iment was performed with a reversed applied magnetic 
field. This seems to imply that the measured dichroic 
signal has a time-reversal odd character (like that of a 
non-reciprocal tensor), contrary to the case of the elec- 
tric quadrupole L 2 — L 2 . In the alternative explanation 
given in Refi^, the change of sign measured in2£, when 
the magnetic field is reversed, is interpreted as due to the 
same electric quadrupole, but in a different twin chrys- 
tallographic domain. In fact, V2O3 has the tendency 
to form big monodomains in its monoclinic phase, in an 
unpredictable way2S, and there are some specific angles 
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at which the measured intensity in one of the other two 
monoclinic domains appears exactly reversed. I shall not 
go into details here, referring to the original papeoS 5 - for 
technical aspects. But I want to underline the following 
two points. First of all, this interpretation is coherent 
with all past failures to detect a magnetoelectric effect in 
V2O5 3334 , and docs not require an ad-hoc mechanism to 
justify a magnetoelectric space group in the ground state 
of V2O3. In second place, such an interpretation can be 
checked by some very simple experiments, that I shall 
describe below. I indeed believe that this is the only way 
to really settle the question. 

There are basically three ways to get some insight into 
the problem. The first, and probably best, is to repeat 
the experiment performed in RefS, but paying particu- 
lar care to check the monoclinic twin after each transition 
to the AFI phase, for example by means of an x-ray- 
scattering equipment: this would give a full geometri- 
cal control of the system. The role of the magnetoelec- 
tric annealing 20-25 , in this sense, should also be analyzed. 
Another possible experimental test to confirm, or reject, 
the interpretation given in Ref^ 5 -, is to measure the lin- 
ear dichroism in the AFI phase in the same geometrical 
conditions as in^fi, ie, with the x-ray direction parallel 
to the corundum axis, but without any magnetoelectric 
annealing, that would be irrelevant to the purpose: in 
this way the El-El origin of the signal can be proved or 
disproved. Moreover, performing an azimuthal scan of 
the dichroism around the corundum axis, it is possible to 
compare directly the signal to the simulations performed 
inS 5 -. Finally, another way to unravel the contradictions 
is to study the dependence of twin formations on the 
applied electric and magnetic fields by means of an x- 
ray-diffraction equipment, in order to determine if their 
presence can influence the phase transition from PI to 
AFI and establish whether the ground-state or some ex- 
cited state properties are investigated by this procedure. 

In the light of all what stated above, and given the 
implications of the non-reciprocal effect, if confirmed, for 
the physics of V2O3, its unambigous experimental deter- 
mination is in my opinion of the utmost importance. 



IV. THE GROUND-STATE IN THE AFI PHASE: 
CORRELATED MAGNETIC MOMENT AND 
NON-LOCAL EXCHANGE. 

In this section I shall focus on two apparently un- 
related facts that are nonetheless founded on a com- 
mon background: i), the reasons why LDA (and related 
LDA+U) are completely unreliable in the description of 
the ground-state wave function of V2O3; and ii), the cor- 
rect symmetry analysis to explain why the magnetic mo- 
ments lie in the glide-plane. The first task is achieved 
through a critical review of all p ie fsilSiiLiSi2£i2&; the sec- 
ond is attained by simply analyzing the local symmetries 
as determined by the molecular nature of the ground- 
state wave function, whose point group does not coincide 



with that at the vanadium site. Both effects are a conse- 
quence of the non-local correlations in the wave-function 
(whose specific form will be precised below) . 

Before starting the theoretical analysis of the ground- 
state wave function, I prefer to introduce the experi- 
mental results at the basis of the molecular model, that 
are probably much more indicativ e) 36 ' 37 . In 1976 J.W. 
Allen 36 wrote that " The magnetic and optical properties 
of all the phases of V2O3 show a loss of V 3+ -ion identity" . 
This statement was based on the optical experimental 
study of (V 2 3 ) x (A\ 2 3 )i- x mixtures (0.01 < x < 0.08). 
The paper continues with: "Optical studies [...] show 
that V 3+ -V 3+ interactions between second and further 
near neighbors are too weak to destroy the V 3+ -ion iden- 
tity. This result suggests the possibility that the basic 
localized-electrons unit in the antifcrromagnetic insulat- 
ing phase may be the nearest neighbor pairs, with the 
electrons delocalized within the pair". I think that it 
can be instructive to compare this statement to the one, 
of opposite tenor, based on LDA(+U) calculations, from 
Elfimov et alm&: "[...] we show that there are other hop- 
ping integrals which are equally important for the band 
shape as the integral for hopping between the partners of 
the pair" . Why do LDA (+U) results differ so much from 
optical data and from all theoretical descriptions^ 5 * 1 ^ 7 - 
? The answer to this question is already contained, im- 
plicitly, in the discussion of section V in Ref>i&, to which 
we refer for a quantitative treatment: here we just try 
to delineate qualitatively the physics behind it, trying to 
add new elements for a better comprehension. As a re- 
sult, we shall be able to define clearly the reasons behind 
the failure of local density approximations to describe the 
physics of V2O3. 

As illustrated in Refii&, a good starting point to eval- 
uate the energetics of the vertical molecule is through an 
orbitally-degenerate Hubbard Hamiltonian. We shall dis- 
cuss a posteriori why it is possible to consider simply the 
vertical pair, that is the contested point in Ref^S, and 
why our results do not depend on the approximations 
that are impicit in the Hubbard model. The parame- 
ters that play the most relevant role are: the Coulomb 
on-site repulsions Ui and Ui-, between electrons on the 
same orbital and on different orbitals, respectively; the 
Hund's coupling constant Jh and the trigonal field split- 
ting A t . The relation U\ = U2 + 2J# holds with very 
good approximation. In V2O3 every V 3+ -ion is octahe- 
drally surrounded by a cage of oxygens and, thus, the 
two valence electrons are forced to a three-fold t 2g con- 
figuration. The effect of the trigonal field splitting is to 
lift the three-fold i2g-degeneracy and to push higher in 
energy the a\ g state, leaving as a ground state the de- 
generate doublet e UtV . The explicit expression of a\ g , 
e u _ v in terms of 3c?-orbital wave-functions can be found, 
for example, ini£. The main contribution to the kinetic 
energy in the molecule comes from the terms describing 
ai g -aig hopping, p ~ -0.82 eV, and e UiV -e UlV hopping, 
fi ~ 0.2 eV— (notice that they have opposite signs). In 
the high trigonal-field-splitting limit the ground state is 
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given, thus, by = \e u e v ; e u e v ) (the two groups of e Uj „ 
labels refer to the electron occupancy of one or the other 
ion in the molecule, e.g., in Fig. V\ or V2). Even for 
the trigonal splitting A t ~ 0.4 eV, deduced from Ref— , 
the ground-state wave function is characterized by a high 
percentage of e ff -doublet composition. Notice that sim- 
ilar results were already obtained in 1978 by CNR, in 
their second paper—, through an unrestricted Hartree- 
Fock + Hubbard U calculation, which is very akin to the 
modern LDA+U. A complete description of the molecu- 
lar behaviour within a big range of variation of parame- 
ters can be found in Ref*i£. Here we focus on the values 
given in Ref— , even if Jh is clearly overestimatedi&^i in 
the light of the atomic calculations of Tanabe and Sug- 
ano {J H = 0.79 eY^-. A t ~ 0.4 eV; J H ^ 0.93 eV; 
U2 — 2.8 eV. Consider now the FM correlated molecular 
wave function 41 -, first derived by Mila et al*^: 



|*o) = —^(\ai g e u ;e v e u ) + \e v e u ;a lg e u ))<3\S M = 2) (3) 

Its energy, in our description, is given by: Eq = 
A t -(p-p) 2 /{U 2 -Jh) ^ -0.16 eV. Notice that the refer- 
ence energy E\ = has been assigned to the ground state 
of the LDA+U procedure ^1 = \e u e v ; e u e v ): here the ab- 
sence of dig electrons gives no trigonal-field contribution 
and Pauli principle prevents any hopping in the ferromag- 
netic vertical molecule. One might wonder why is Eq so 
lower, having used the same parameters as obtained by 
the local density calculations. To this aim, consider the 
wave function | ^2) = |e„e„; ax g e u ). The associated en- 
ergy is given by: E 2 = A t - (p 2 +p 2 )/{U 2 - J H ) ^ +0.02 
eV, higher than both Eq and E\. There are 180 meV of 
difference 4 ^ between the "mean-field" state \^ 2 ) and the 
"correlated" state l^o)* coming from the hopping inter- 
ference term 2pp/{U 2 — Jh)- The origin of this energy is 
due to the correlated structure of |^o) that makes pos- 
sible for an e„-electron on Vi to hop on V 2 and simul- 
taneously for an a\ g electron on V2 to hop on Vi (and 
vice- versa), without changing state. This form of su- 
perexchange interaction is basically a non-local electronic 
correlation and, as such, it cannot be taken into account 
in the framework of local density calculations. These sim- 
ple considerations show why LDA+U results fail to catch 
the ground-state wave function: as non-local correlations 
are neglected, the interference energy —2pp/ (U 2 — Jh) is 
completely missed (~180/2=90 meV per ion !) and the 
purely e UyV state is found as ground state. 

In this picture, if we even consider that there are three 
in-plane second nearest neighbors and introduce further 
nearest-neighbor hoppings, as done in Ref— , still the 
non-local correlation energy of the vertical molecule rep- 
resents a relevant (if not the main) contribution. Notice 
that the corundum geometry is such that further-nearest- 
neighbors hopping integrals cannot give rise to any other 
correlated hopping of relevant strenght: for this reason 
all inter-molecular interactions can be treated in mean 
field (or local density) approximation. The energy gain 



due to the correlated molecule is so big that it can be 
hardly conceived as an artefact of the Hubbard Hamil- 
tonian, or that the corresponding state dissolves when 
further nearest neighbor hopping-integrals are taken into 
account, as Ref— seems to suggest, provided, of course, 
all local and non-local forms of energy are correctly con- 
sidered. Moreover, it is stable against rather large varia- 
tions of A t , U 2 and Jh around the chosen values— When 
this correlation energy is simply neglected, as in local 
density approximations, the ground-state wave function 
is mainly determined by the trigonal splitting A t , that 
pushes higher the one-particle energy level of the a\ g - 
electrons, thus advantaging e UjtI . In my opinion, this 
is a signature of the failure of LDA+U to describe the 
ground state of V2O3. Indeed, the results of LDA+U 
calculations^ are lacking for many aspects: the spin ex- 
changes along the different bonds are strongly underesti- 
mated when compared to the experimental results^, as 
recognized by the authors themselves. The orbital occu- 
pancy of mainly e u e v kind does not reproduce the find- 
ings of Park et a/— interpreted with a substantial a\ g - 
occupancy about 20-25%. Notice that this is in keeping 
with the experimental evidence from Rubinstein 44 that 
found no detectable anisotropy in the Knight shift of the 
NMR spectrum, no observable quadrupole splitting in 
the NMR spectrum, and zero anysotropy in single-crystal 
susceptibility measurements for various cristalline orien- 
tations with respect to the magnetic field. All these find- 
ings apparently point to a negligeable trigonal field split- 
ting, ie, an almost cubic crystal field. Only the molecular 
picture, through the interference effect, is able to recon- 
cile these latter experiments with the sizeable value of 
the trigonal distortion, because of the average orbital oc- 
cupancy ai g = 25%, e MjtI = 37.5%, close to the statistical 
ratio (33%) of the cubic field. In the light of the previous 
discussion, the main statement of Elfimov et aL 18 looses 
its strenght: these authors performed a tight-binding fit- 
ting of their LDA and LDA+U band structure, and con- 
cluded that a proper description requires second, third 
and forth neighbor hoppings. Yet, the evident drawback 
of LDA+U correlated description makes such a conclu- 
sion very weak. Moreover, they focused on the a\ g band 
only, whereas the three bands should strongly interfere. 

My conclusion in this respect is that V2O3 is 
a very peculiar strongly correlated electron system, 
where ordinary ab-initio codes may loose their original 
strenghtpi&^Si— The case of LDA+U is emblematical: 
this approach was initially conceived to describe the in- 
sulating magnetic materials whose energy gap was not 
determined by the Stoner parameter, but by the Coulomb 
repulsion £/— . Even though it works correctly in these 
cases^, it is nonetheless based on a local exchange and 
correlation potential, and, as such, it cannot treat prop- 
erly a system like V2O3, where non-local correlations 
amount to 90 meV per ion. The same arguments apply 
to the recent LDA+DMFT 4 &, as clear from the following 
statement, from Keller et al&. "[...] this small trigonal 
splitting strongly determines the orbital ground state of 
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the V-ion obtained from LDA+DMFT calculations [...]". 
This method, in spite of the progress coming from the 
DMFT, can handle just the local dynamics of the sys- 
tem, while non-local correlations are still neglected. This 
is the reason why the orbital nature of the ground state 
seems to be determined by the trigonal field splitting. 

It may be true that the physical idea of the " weakly" 
interacting gas of vertical molecules^ could be oversim- 
plified and that a different picture for the band structure 
may be required, as in the original idea of Elfimov et ali&: 
yet, to this aim, an extension of the correlation potential 
V c in a non-local direction is, in my opinion, essential to 
describe the physics of V2O3. All past experimental re- 
sults, like Allen's, must be explained by any theory, and 
not simply dismissed, as seems to be the case of LDA- 
based approaches. In order to get reliable results, it is 
necessary to go back to the original form of Kohn-Sham 
equations^ and find out a good approximation for the 
correlation potential V c : some hints in this sense could 
be taken from the form of the non-local molecular wave 
function (J3J. In this way new informations about V2O3 
band structure can be extracted, that do not neglect the 
"correlated hopping" energy, and, probably, can also ex- 
plain, among the others, Allen's results^. 

All the previous analysis was based on the non- 
local nature of superexchange correlations, only, and did 
not take into account the spin-orbit interaction, that, 
nonetheless, was fundamental to explain the results of 
Notice that, even if very important to de- 
scribe the RXS experiments and to remove the orbital 
degeneracy, this interaction is about one order of magni- 
tude lower than the correlated hopping (~ 25 meV) and, 
thus, it can be introduced adiabatically in the ground- 
state manifold of | ^o) • I n this way, its inclusion, though 
it removes the orbital degeneracy of l^o) 1 — , leaves unal- 
tered all the previous conclusions. The ground state is 
modified in: 



|*o) = -j={\a lg e + -e+e ) + \e+e ; a lg e + )) <g> \S Mz = ±2) 

(4) 

where e + = — "^75 ( e u + ie«) and e = -^75 (e« — ie v ). Here 
the spin part of the ground state has been written ex- 
plicitly. Such a state is characterized by a global spin 
Sm = 2, with Smz = ±2. As discussed previously, a long 
debate arose in the literature about a possible component 
of the spin of V2O3 out of the glide plane, especially for 
its consequences on the magnetic space group. From a 
theoretical point of view, on the basis of usual crystal 
field models, such an out-of-plane component could even 
be expected: in fact, vanadium ions, in their AFI phase, 
are characterized by the absence of any point symme- 
try, thus allowing the local magnetic moment to follow, 
in principle, any direction. Yet, the molecular nature 
of the ground state plays a fundamental role also in the 



description of the correct symmetry. In fact, as any phys- 
ical variable is determined through its expectation value 
on the ground-state wave function, and as the electronic 
wave function is " molecular" , then the magnetic moment 
must be invariant with respect to all the symmetry op- 
erations of the molecular point group and not of the ion 
point group. In other words, we identify the molecular 
spin Sm = SVi + SV2 = 2 as the physical variable de- 
tected by neutron measurements, even if it is customary 
to speak about the magnetic moment per ion, Sy = 1, 
in an average sense. This interpretation parallels and ex- 
tends to the four-electron case of the molecule the Sy=l 
atomic picture: here, of the whole space spanned by the 
two S = 1 /2-electrons on the vanadium ion, the Sy=0 
subspace is projected out, and only the spin-1 is consid- 
ered for the ground-state properties. Similarly, in the 
molecule, only the Sm = 2 subspace has a physical re- 
ality, and all the other spin-states (Sm = 0, 1) are pro- 
jected out. 

Based on these premises, we are now ready to repeat 
the previous symmetry considerations from a different 
viewpoint. In fact, the electronic ground-state of the 
molecule now has a local symmetry element, contrary to 
the case of the single ion: this symmetry element is the 
C2b-axis shown in Fig. ^ orthogonal to the molecular 
axis, and running through the midpoint of the molecule. 
It follows that the molecular magnetic moment must lie 
either along the C2-axis, or perpendicular to it (this latter 
case arises when the time-reversal operation is associated 
to the C2-axis). No "mixed" configurations are possible. 
Neutron experiments allow us to choose: the magnetic 
moment lies orthogonal to the C2j>-axis and, thus, the 
correct local symmetry for the molecule is TC 2 . The 
same conclusion is reached by considering Tanaka's wave 
function In fact, under C 2 symmetry, a\ g — > ai g , 
e + — > e + and e~ — > — e~ (seeii, pag. 3), and time- 
reversal symmetry changes sign to the imaginary unit. 
Thus, (%\L M \%) = (%\(fC 2 )^L M fC 2 \%)^nd the 
molecular angular moment is invariant under TC 2 sym- 
metry, so that it must lie on the glide plane and force 
the spin moment in the same direction. Notice, in this 
respect, that when Tanaka speaks of a possible out-of- 
plane magnetic component he refers to the single-ion mo- 
ment: when combined throughout the unit cell, all the 
ions belonging to the same vertical molecule have oppo- 
site out-of-plane moments and thus the global molecular 
spin is always in the glide-plane, as the symmetry con- 
straint requires. 

In conclusion, I have shown that the non-local, molecu- 
lar, wave function can explain different physical phenom- 
ena, apparently unrelated and inexplicable in the frame- 
work of an average, ionic picture. 

I would like to thank Stephen W. Lovesey who strongly 
supported and stimulated the present work, and Calogero 
R. Natoli who shared with me his knowledge on V2O3. 
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